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derive general properties of the critical curves and caustics analytically. The 

r> ' multiple-image cross section is a sensitive function of the inclination angle of 

the disk relative to the observer. We compute the inclination averaged cross 

section for several sets of lensing parameters. For realistic disk mass and size 

parameters, we find that the cross section for multiple imaging is increased 

by only a modest factor and no dramatic increase in the optical depth for 

strong lensing of QSOs would be expected. However, the cross section for high 

magnifications is significantly increased due to the inclusion of a disk, especially 

for nearly edge-on configurations; due to the strong observational selection 
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effects favoring high magnifications, there might be significant consequences for 
lensing statistics. 
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1. Introduction 

Early theoretical studies of gravitational lensing of distant quasars by foreground 
galaxies approximated their mass distributions as spherical and usually also as singular and 
isothermal QGott fc Gunn 1974| ; p?urner 1980| ; Turner, Ostriker & Gott 1984). Later studies 



explored the relaxation of these simplifying approximations and, in particular, considered 
elliptical mass distributions and potentials {e.g., piandford fc Kochanek 1987| , [Kochanek 



fc Blandford 1987| , [Blandford fc Narayan 199"^ ), but it has been shown that the standard 



spherical approximations are adequate for the purpose of many statistical calculations {e.g., 
Fukugita fc Turner 1991| , |Maoz fc Rix 1993| , [Kochanek 1993|) . However, with rare exception 



( Pstriker fc Vietri 1990| ), the gravitational lensing effects of the defining component of 
spiral galaxies, their very thin disks, has been ignored in discussions of lensing statistics. 
There are two reasons. First, within the limitations of the spherical galaxy approximation, 
it was early shown (Turner et al. 1984) and later repeatedly confirmed by more detailed 
studies ( [Fukugita fc Turner 199l| , [Kochanek 1991[) that spiral galaxies contribute only 



a small fraction, of order 10%, of the total lensing cross section of the known galaxy 
population, with ellipticals and SO's being the dominant contributors. Second, the disk 
stability arguments originated by Ostriker & Peebles (1973) and the many other lines of 
evidence which indicate the presence of massive dark halos suggest that even spiral galaxy 
mass distributions are dominated by a roughly spherical component. 

While these considerations remain valid, there now seem to be some reasons to examine 
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the lensing effects of thin disk mass distributions more carefully. For one, samples of 
galaxy lensed quasars are becoming rapidly larger and better controlled. Thus, smaller 
and more subtle features of lens statistics are becoming observationally accessible and 
therefore theoretically interesting. For another, in the usual (and very accurate) thin lens 
approximation, it is not the three dimensional mass distribution of the lensing object that 
matters but rather the projected two dimensional one, particularly those regions which 
exceed the lensing critical surface density ([Turner et al. 19841) . From this two dimensional 



point of view, a very thin disk may be much more important when seen nearly edge-on 
than it is in three dimensions. Finally, although the situation is far from unambiguous, 
observations suggest that lensing galaxies may have more dust ( [Lawrence et al. 1995| , 



Malhotra, Rhoads & Turner 1997) and larger quadrupole mements (Keeton, Kochanek 
& Seljak 1996) than would have been naively expected if early type galaxies entirely 
dominated the quasar-galaxy lensing rate. 

Therefore, in this paper we present a preliminary exploration of the lensing properties 
of spiral galaxies, which we approximate as finite, uniform surface density disks of zero 
thickness embedded in singular isothermal spheres (SIS). We focus on the general caustic 
and critical line properties of such lensing objects and also study the lensing cross section 
enhancement over the SIS only approximation. We explicitly carry the important variable 
of disk inclination angle throughout the calculations. In the present paper, however, we do 
not present the fully detailed calculations required for direct comparison to observations; 
in particular, we do not here consider line-of-sight integrations, integration over realistic 
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distributions of galaxy properties nor the effects of amplification (magnification) bias on 
observed samples. Rather, we present only an initial qualitative and semi-quantitative 
investigation of lensing by disk galaxy mass distributions. 

The basic equations and model are presented in section 2. Section 3 discusses our 
technique for finding caustics, critical curves and multiple imaging configurations. Section 
4 presents the general lensing properties of our model of disk galaxies in the language of 
caustics and critical curves. Section 5 considers the effects of inclined disks on lensing cross 
sections and the implied modification of quasar lensing rates. Finally, section 6 contains a 
summary of our results and a discussion of recent related work. We follow the conventions 
and notation of Schneider, Ehlers & Falco (1992). 

2. Light deflection due to a spiral galaxy 

We model spiral galaxies as uniform and infinitely thin disks embedded in singular 
isothermal spheres. Let us choose a length scale appropriate for a singular isothermal sphere 
(SIS), 

«o^4.(-H)^i^. ,U 

where vsis is the line-of-sight velocity dispersion of the SIS. Let $, and rj be the physical 
position vectors of the image (in the lens plane) and the source (in the source plane) 
respectively, then the dimensionless image and source positions are 

x = f; y = ^. (2) 



-6- 

where tjq = ^qDs/ Dd, Dg and D^ are our distances to the source and lens respectively. The 
lens equation becomes 

y = x-a(x), (3) 

The scaled deflection angle 



- /dV«:(x') A^' 
n J X — x' 



a(x) = £/dV«:(x')^:^^, (4) 



where k(x) = S(^ox)/Scrit, with T^crit = c^Ds/^AnGDdDds)- Dds is the distance between 
the lens and the source. We use afline distances, Dd = cHQ^\{zd), Dg = cHq^\(zs), 
Dds = cHo^\{zd, Zs) = cHq'^{1 + Zd)[\{zs) - \{zd)], where 



A(;.) = r dw , ^ (5) 



For Qo = 1, X{z) = (2/5) 



1 - fl + z)-^/^' 



For a SIS, its dimensionless surface density and scaled deflection angle are ( [Schneider 
et al. 19921 ) 



'^"^"(x) = ^; a^^^(x) = ^. (6) 

2|x| |x| 

The projection of a circular disk of radius rdisk in the plane perpendicular to the line 
of sight is an ellipse with semi-major axes rdisk and rdisk cos 9, where 9 is the angle between 
the normal vector of the inclined disk plane and the line of sight. Note that < 9 < 7r/2; 
6' = is the "face-on" disk case; 9 = tt/2 is the "edge-on" disk case. 

For a uniform disk, its projected surface density is T,disk = Mdisk/{'^i^disk'^'^^^)j i^^ 
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projected dimensionless surface density is 



disk 
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K 



disk 






risk COS 9 
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cos 6* V Uisk . 



^GMd^skPdPds 



D. 



D 



ds 



1000 Mpc/ Ds 



(7) 



For simplicity in notation, we define 



q = /t**''cos6'. 



q is the face- on surface mass density of the disk. Note that g < 1 for reahstic choices of the 
parameters. 



The dimensionless radius of the disk is 



R 



r-disk 10 / r-d, 



isk 



'I50km/s\ /1000Mpc\ D 



^0 TT \ 10 kpc 



vsis 



V Dd 



D. 



(9) 



ds 



Let us choose the x coordinates such that the contour of the projected disk is given by 



r) [rcosO 



1. 
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Then the total scaled deflection angle a = cx'^^^ + a , with a given by 



a 



disk 



duj In 



an 



, disk 



Rq 

Rq 

2ticos9 7-1 



(Vl-tu2 + xi/Rf + (a; cos 9 - X2/Rf 



(VI -^^ - xi/Rf + {uj cos 9 - X2/Rf 



duj In 



{U - Xi/Rf + (VI - UJ'^ cos 9 + X2/Rf 

_(CJ - Xi/RY + (VI - cc^2 cos^ - X2/RY_ 



:iii 



Note that 



ttffc(cos^^O) = ^ C duo 



TX 



V^ 



ijj^ 
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To study image multiplicity, we will also need the derivatives of the deflection angles: 
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We have defined 

/l(cc') = 
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cucos^ 



iU \ + I — — Vl — CJ^ COS ^ 
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3. Condition for multiple images 

The magnification factor of the source is given by /i(x) = 1/det A(x), where the 
Jacobian matrix A(x) is defined as 



A(x) 



detA 



dy 



A. 



1- 



dai 
dxi 



1 



(9a;j ' 



(15) 



9X9 



9q;i (9q;2 
9x2 9x1 



It has been shown that an isolated transparent lens can produce multiple images 
if, and only if, there is a point x with det A(x) < 0. If at xq, det A(xo) < (negative 
parity), a source at yo = y(xo) has at least two additional images of positive parity. The 
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multiple-iinage cross-section is simply the bounded area in the source plane in which for 
each source position y, there exists an image position x where det ^(x) < 0. 

Formally, critical curves are given by det A{x.) = 0, the corresponding source positions 
are the caustics. However, the cross-section for multiple images is not always bounded 



by a caustic. In the case of SIS only, det A = 1 — 1/x, where x = Jxl + xl- det A < 
gives < X < 1. The lens equation gives ns y = |x — 1|. Hence y < I for multiple images; 
the multiple image cross-section is bounded by y = 1, which corresponds to x = (where 
det A = — cxd), while the only critical curve is at x = 1 corresponding to y = 0. 

Next let us consider the case of SIS plus a face-on disk, ^ = 0, k'^**^ = q. It's 
straightforward to integrate Eqs. fpUD to find ( Schneider et al. 1992|) 



a 



x/x + K X, X < R] 

x/x + k'^^'^'^R'^ x/x^, X > R. 



(16) 



Hence we have 



det A 



1 - k'^''^^ (^1 - /t*^'^ - 1/x) , X <R; 

'l + k'^'^^R^/x'^) (l - 1/x - k'^''^R^/q 
It is easy to see that det A is discontinuous at x = R. 



(17) 



'X 



X > R. 



For k'^**'^ < 1, detA does not change sign at x = i?. If i? > 1/(1 — k*'^'^), the only 
critical curve is given by x = 1/(1 — k'^*'''^) < R, and it maps to y = 0; if i? < 1/(1 — k'^**'^), 
the only critical curve is given by x = f 1 + Vl + Ak'^^^^R'^j /2 > i?, and it also maps to 
y = 0. The critical curves therefore have no relevance to image multiplicity. On the other 
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hand, det A{x) < leads to y < 1, with y = 1 given by x = 0, just as in the SIS only case. 
SIS with a face-on disk of k'^*'^^ < 1 has the same multiple-image cross section as SIS only. 

If k''-'^^'' > 1, the true critical curve is given by a; = (1 + Vl + Ak'^^^^B?)/2 > R, which 
again maps to y = 0. However, det A{x < R) > 0, det A{x — » R^) < 0, hence det A changes 
sign at X = i?, which is effectively a "critical curve" ; the corresponding "caustic" is given 
by y = 1 + R[k'^^'^^ — 1) > 1. SIS with a face-on disk of k'^*'^^ > 1 has a larger multiple- image 
cross section than SIS only. 



4. Image multiplicity for SIS plus inclined disk 

For SIS plus uniform disk, a.'^'^^^^x = 0) = 0, hence a; = corresponds to source position 
y = 1, just as in the SIS only case. When det A{x ^ 0) < 0, y < 1 always gives multiple 
images; the multiple-image cross-section increases relative to the SIS only case if any caustic 
curves lie outside the y = 1 circle. 

The value of det A at x = is a useful indicator of the general properties of the critical 
curves. For uniform disk plus SIS, we find [using Eqs.([T3l) and ([T5|) ] 

2g , 1 / 2g \2 C,{e) xl , C2{9)xl 



det A(x ^ 0) = 1 ^ + 



cos^ ' cos^ Vl + cos^y x^ x^ 

where 

l + cosfc' costal + costal 
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Let us define 

S = Cixl + C2XI. (20) 

S has the same sign as det A{x —>■ 0). We have two critical angles: 

Ci{9) = 0: 9 = 9i = aTCCOs{2q - 1); 



C2{9) = : 9 = 92 = arccos [ V^+^g M _ (2I) 



Note that 9i is not defined for g < 1/2. For g > 1/2, ^1 > ^2- 

Note that det A is discontinuous for (x5,a;2) on the ellipse {x\/ RY + {x2/Rcos9)'^ = 1 
[see Eqs. QI^) ], i.e., on the edge of the disk, just as in the face-on disk case. If 
detA{x -^ 0) = +00, det y4 changes sign at (x^,X2), which describes effectively a "critical 
curve" , and the corresponding source positions give us a "caustic" . 

First we consider q < 1/2. Here we always have Ci < 0. We have two different cases: 
(1) ^ < ^2, C*2 < 0, S* < 0, detA{x — » 0) = —00, we have only one true critical curve and 



it encloses x = 0; (2) 6* > 6*2, C2 > 0, S > [detA{x -^ 0) = +00] if \x2\ > J\Ci\/C2 \xi 



S < [detA{x — > 0) = —00] if \x2\ < J\Ci\/C2 \xi\, here we have one true critical curve 
enclosing x = 0, plus an additional "critical curve" which is hourglass shaped and passes 
through X = 0, and is closed off by {x'i,X2) defined above. Figs. 1-3 show the critical curves 
(a) and caustics (b) for q = 0.2 (^2 = 72.168°), for 9 = 70°, 75°, and 85° respectively. The 
unit radius circle in (a) indicates the critical curve for the SIS only case (its corresponding 
caustic shrinks to the point y = in the source plane); the unit radius circle in (b) indicates 
the multiple-image cross-section for the SIS only case. 
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For q > 1/2, we have three different cases: (1) 9 < 02, Ci < 0, C2 < 0, S < 0, 
det A{x -^ 0) = —00, we have only one true critical curve and it encloses x = 0; (2) 



^2 < ^ < ^1, Ci < 0, C2 > 0, 5 > [detA{x ^ 0) = +00] if \x2\ > \\Ci\/C2 |xi|, 5 < 



[det 74(a; ^ 0) = —00] if 1x2! < J\Ci\/C2 \xi\, here we have one true critical curve enclosing 
X = 0, plus an additional "critical curve" which is hourglass shaped and passes through 
X = 0, and is closed off by (x^jXg) defined above; (3) 9 > 9i, Ci > 0, C2 > 0, S > 0, 
det A{x ^ 0) = +CX3, we have two critical curves enclosing x = 0, one of them is the true 
critical curve, the other is given by {xl,X2) defined above. Figs. 4-6 show the critical curves 
(a) and caustics (b) for q = 0.6 (^i = 78.463°, 62 = 45.238°), for 9 = 45°, 75°, and 85° 
respectively. The unit radius circles in the figures are the same as already indicated. 

Because det A{x) is modified from Eq.(r3) as we move away from a; = 0, the hourglass 



shaped "critical curves" are given by the intersection of \x2\ ~ \/\Ci\/C2 \xi\ and the ellipse 
{xI/RY + {x^l Rcos9)'^ = 1; this explains the curved sides of the hourglass shape. Due 
to numerical noise, some discrete points from the {xl,X2) ellipse and their corresponding 
source positions appear in Figs. 1-6 as discrete points apart from the critical curves and 
caustics. 

It is interesting to note that our study of the face-on disk case (see the previous section) 
suggests that the critical angle for the disk to contribute to image multiplicity is given by 
^disfc _ q,^ cos^ = 1, i.e., 9 = 9q = arccos(g). For q < 1/2, only 6*2 is defined, and 9q > ^2! for 
q > 1/2, ^1 > ^0 > ^2- This indicates that an inclined disk is more efficient than a face-on 
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disk with the same surface mass density. 



5. Increased multiple-image cross section due to the inclined disk 

When there is only one critical curve (the true critical curve) in the x plane, the 
multiple image cross-section is given by the area enclosed by y = 1 plus the areas enclosed 
by the caustic which lies outside the y = 1 circle, which are two small sharp-angular areas 
lying along the t/i axis; the caustic is completely inside the y = 1 circle for sufficiently 
small 9. When we have two branches of critical curves (one of them is a true critical curve, 
the other corresponds to discontinuity and change of sign in det A) in the x plane, the 
multiple image cross-section is given by the area enclosed by y = 1 plus the areas enclosed 
by the caustics which lie outside the y = 1 circle, which now consist of two angular areas 
lying along the yi axis (corresponding to the true critical curve, at largest possible x), 
plus two round areas lying along the 1/2 axis (corresponding to the "critical curve" due to 
discontinuity and change of sign in det A). 

Assuming uniform distribution of the inclination angle 6 in solid angle, the average 

multiple-image cross section is 

Jd9sm0aie) 
J dO sm 

where a{6) is multiple-image cross section for inclination angle 6. 

Figs. 7-9 show the ratio of the multiple-image cross-sections for SIS plus inclined 
uniform disk and for SIS only, a{9)/o"^^^ , as a function of the inclination angle 6, for three 
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groups (nine sets) of choices of (g, R) . Table 1 lists the parameters and the corresponding 
average enhancement factor in multiple-image cross section ag/a^^^. 

Table 1: Average enhancement factor in multiple- image cross section 





R= 1.5 


R = 3 


R=6 


q=0.2 


1.086 


1.160 


2.041 


q=0.4 


1.296 


1.547 


2.742 


q=0.6 


1.718 


2.493 


4.509 



Let us write 



qR 



1.9 flOkpc 



M, 



disk 



150km/s^ 



(23) 



TT V ^disk J \10^^ Mq J \ VsiS J 

For a model spiral galaxy with given [vsis, M^isk, fdisk), q Ris constant. Inspection of Table 
1 shows that the modification of the galaxy's inclination averaged cross section is quite 
small, less than say 50% (there are other uncertainties in lens statistics calculations at least 
this large), if gi? is less than about unity. From equation (23) we then see that significant 
modification of the multiple image lensing cross section will only occur for objects with 
uncharacteristically small and massive disks and/or those with minimal spherical (halo) 
components. Few, if any, real galaxies may satisfy such conditions. 



The differential probability of a strong-lensing (multiple-images) event is ( p?urner et al. 



19841) 



dr = riLizd) (7e{zd\zs) 



cdt 

dzd 



dzd 



(24) 
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where ni^z^) is the number density of lenses at lens redshift z^, a0{zd\zs) is the inclination- 
angle averaged cross-section for multiple images given source redshift Zg- From the above 
discussion, we would expect the total contribution of spiral galaxies to the QSO strong 
lensing optical depths to increase only slightly due to the effects of their thin disks. The 
optical depths would then continue to be dominated by early type galaxies. 

However, the inclusion of an inclined disk breaks the circular symmetry due to the 
SIS, the true critical curve [where detyl(x) = 0] no longer maps to a point {y = 0) as in the 
SIS only case, but maps to a caustic which encloses an area comparable to the SIS only 
multiple-image cross-section for 6 > 62 [see Eq.(pT|)]. Since the magnification is infinite (for 
a point source) on the caustic, and decreases smoothly away from it, the cross section for 
high magnifications is significantly increased due to the inclusion of a disk. Quantitative 
investigation of this effect will require extensive numerical calculations which are outside 
the scope of this paper, but this effect may well be important. Observational selection 
effects favoring inclusion of high amplification (magnification) lensing configurations in flux 
limited samples can lead to major distortions of intrinsic distributions in real samples; see 
Ostriker & Vietri (1990), for example. 



6. Discussion and conclusion 

We have examined strong (multiple image) gravitational lensing by spiral galaxies, 
modeled as infinitely thin uniform disks embedded in singular isothermal spheres. We 
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have derived general properties of the critical curves and caustics analytically. The 
multiple-image cross section is a sensitive function of the inclination angle of the disk 
relative to the observer. We have therefore computed the inclination averaged cross section 
for several sets of lensing parameters. 

We find that the optical depth for multiply imaged QSOs should only increase by a 
factor of a few at most and by less than 50% in nearly all realistic cases; the inclusion of a 
disk is therefore expected to have no significant effect on the contribution of spiral galaxies 
to the total optical depth for multiply imaged QSOs. On the other hand, the cross section 
for high magnifications is significantly increased due to the inclusion of a disk. The increase 
in the number of lensed QSOs with high magnifications could have substantial effects on 
observed lens samples. 

While completing this work, we became aware of recent preprints by Loeb (1997) and 
by Mailer, Flores, and Primack (1997). The former is primarily concerned with the possible 
connectin between lensing by spiral galaxies and high column density HI absorption systems 
seen in QSO spectra but also briefly mentions the potential role of disk inclination in spiral 
galaxy lensing statistics. The subject of our paper overlaps partially with that of Mailer et 
ai, and some of our results are qualitatively similar; however, our paper is mainly analytical 
(which provides insight into the mathematics of the critical curves and caustics) while theirs 
is numerical. 

We thank A. Loeb for useful discussions and gratefully acknowledge support from NSF 
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Fig. 1.— q = 0.2, i? = 3, ^ = 70° < ^2 = 72.168°. (a) Critical curve. The unit radius circle 
indicates the critical curve for the SIS only case (its maps to y = in the source plane), (b) 
Caustic. The unit radius circle indicates the multiple-image cross-section for the SIS only 



case. 
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Fig. 2. — q = 0.2, R = 3, 9 = 75° > 62 = 72.168°. The unit radius circles are the same as in 
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Fig.l. (a) Critical curve, (b) Caustic. 
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Fig. 3. — q = 0.2, R = 3, 6 = 85° > 62 = 72.168°. The unit radius circles are the same as in 
Fig.l. (a) Critical curve, (b) Caustic. 



C\2 
X 








xl 



yi 



Fig. 4. — q = 0.6, R = 3, 6 = 45° < 62 = 45.238°. The unit radius circles are the same as in 
Fig.l. (a) Critical curve, (b) Caustic. 
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Fig. 5.— q = 0.6, R = 3, 02 = 45.238° < ^ = 75° < ^i = 78.463°. The unit radius circles are 
the same as in Fig.l. (a) Critical curve, (b) Caustic. 
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Fig. 6. — q = 0.6, R = 3, 9 = 85° > 6i = 78.463°. The unit radius circles are the same as in 
Fig.l. (a) Critical curve, (b) Caustic. 



-23- 



I— I 

m 

b 



Qi 
^ 



8 



6 - 



4 - 



2 - 







T I I r 



T 1 I r 



T I r 



q=0.2; R=1.5, 3, 6 (down to up) 



60 




J I I L 



J I I L 



J I I L 



70 



80 



90 



^(degree) 



Fig. 7. — The ratio of the multiple-image cross-sections for SIS plus inclined uniform disk 
and for SIS only, a{6)/a^^^ , as a function of the inclination angle 6, for q = 0.2, R = 1.5, 3, 



6. 
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Fig. 8. — Same as Fig. 7, for q = 0.4, R = 1.5, 3, 6. 
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Fig. 9. — Same as Fig. 7, for q = 0.6, R = 1.5, 3, 6. 



